A chemically-driven finite-wavelength instability in 
ion-irradiated compound semiconductors 

Scott A. Norris 
March 6, 2013 

Southern Methodist University 
Department of Mathematics 
3200 Dyer Street 
Dallas, TX, 75275-0156 

Abstract 

In recent years, observations of highly-ordered, hexagonal arrays of self-organized nanostructures on 
solids irradiated by ion beams at normal incidence have excited interest in this phenomenon as a potential 
route to high-throughput, low-cost manufacture of nanoscale devices. Experimental and theoretical work is 
converging on the hypothesis that these ordered structures only appear when irradiating binary, or impurity- 
laced monatomic targets. An important question now is to identify the fundamental mechanism behind 
structure formation, and recently many studies on this topic have appeared. 

Two theories have recently been separately proposed. One well-studied analytical model predicts or- 
dered patterns driven fundamentally by a morphological instability. However, such patterns are also ob- 
served in numerical simulations for which the driving mechanism is a compositional instability. To identify 
testable differences between the two theories, we propose here a generalized framework containing pa- 
rameter regimes associated with both previously-proposed mechanisms. Analyzing the full model, we find 
that the compositional mechanism also contains all of the necessary ingredients for ordered patterns, and 
is clearly differentiated from the morphological mechanism by differences in parameter values and predic- 
tions on experiment. This suggests the utility of specific atomistic and experimental measurements to guide 
continued theoretical progress. 



1 Introduction 



Among the many nanoscale patterns that have been observed on ion-irradiated surfaces HJ EL the discovery 
of hexagonal arrays of high-aspect ratio dots on Ar + -irradiated GaSb [3 1 has sparked a flurry of experimental 
and theoretical study into spontaneous pattern formation as a potential route to "bottom-up" fabrication of 
nanoscale devices. The mathematical study of such patterns is a mature field known as bifurcation theory, and 
the ingredients necessary to produce long-range hexagonal order are well known H • First, one needs 
isotropic physics, which is readily achieved by irradiating the sample at normal incidence, or by rotating 
it during off-normal bombardment. Second, one needs a destabilizing mechanism that produces a finite 
wavelength bifurcation - a transition in the linearized equations from stability to instability that occurs first 
for perturbations of a finite wavelength. Unfortunately, none of the early models of ion irradiation predicted 
such a finite-wavelength bifurcation 0[7]|, suggesting a need for the incorporation of additional physical 
effects. 

Two candidates that initially seemed promising were a "damping" term proposed to stand for long- 
range atomic redeposition [8] and biaxial compressive stress injected by the ion beam (see |j6l [TJ for a 
comprehensive discussion). However, recently both mechanisms have been ruled out as sources of instability 
J9j HOB . In parallel with these theoretical struggles, persistent difficulty obtaining consistent results between 
labs has led to a re- examination of the experimental results. It was found that dots, originally also observed 
on irradiated pure silicon Hill 11211 . disappeared from that system when impurities and geometric artifacts 
were carefully removed IfTl 1131 H4ll , and conversely, that the systematic addition of impurities to a clean 
system causes dots to appear H151 1161 1171 H8H . This has led to the strengthening conclusion that two- 
component materials are the necessary ingredient to produce ordered quantum dots H19II . 

This growing consensus has resulted in a number of theoretical treatments of irradiated binary materials, 
describing both the morphology and the concentration evolution. Built on the pioneering work Shenoy Chan, 
and Chason H20L Bradley and Shipman (BS) have recently introduced such a theory for binary materials 
exhibiting the first physically-grounded finite-wavelength bifurcation 11211 122 II , which will be discussed in 
greater detail below. This result has sparked great interest, inspiring additional studies on normal-incidence 
irradiation of two-component materials [23], including a variant where the second material is not originally 
present, but instead co-deposited simultaneously with the ion irradiation H24H , and a similar framework 
for the study of ion-assisted deposition, where both materials are deposited along with the ions, and the 
deposition rate exceeds the sputtering rate 112511 . 

In addition to dealing in some way with two materials, all of the papers cited in the previous paragraph 
share an important common trait: they assume the presence of ordinary, Fickian diffusion, in which atoms 
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of each species diffuse from regions of high concentration to regions of low concentration. However, this 
assumption is not obviously valid: preferential sputtering of one material can lead to residual concentrations 
that are not energetically favorable. In this case spinodal decomposition is generically expected if the atoms 
are sufficiently mobile, and although temperatures of ion-irradiated targets are too low to allow thermally- 
driven segregation, it is known that ion irradiation enhances atomic mobility by several orders of magnitude 
12611 . Additional motivation to consider this effect comes from the significant restrictions on (unknown) 
parameter values required by the BS theory for the production of a finite-wavelength instability. Finally, Le 
Roy and co-workers have recently observed intriguing Ga-enriched 'caps' atop nanostructures on irradiated 
GaSB H27L which were then observed with qualitative agreement in a numerical simulation of the phase 
separation process I128H . 

Unfortunately, comparison of the latter work on phase separation to the BS model is hampered by the 
fact that the former study was performed numerically, rather than analytically, with no analysis of the linear 
instability leading to pattern formation. Hence, in the present work, we propose a generalized framework 
in which both linear instabilities can be analyzed, and in which the associated implications for experiment 
can be compared. Our central result is that spinodal decomposition allows a second finite-wavelength bifur- 
cation, distinct from the bifurcation of Bradley and Shipman, that places fewer restrictions on the unknown 
parameters and can be clearly distinguished from it by a 180-degree shift in the relative phases of the height 
field and the concentration field. This result suggests that both atomistic studies of specific parameter val- 
ues, and experimental measurement of the composition profile at the early stages of dot formation, would 
be equally valuable guides to further theoretical progress. 



Although the ultimate goal is to understand the formation and properties of high aspect ratio dots at late 
experimental times, in this work we will focus exclusively on the linear problem associated with small am- 
plitudes and very early times in the pattern-forming process, for two reasons. First, a nonlinear analysis can 
only be as accurate as the linear instability out of which it grows, and experimental confirmation of linear 
models is still being sought. Second, once a long-wavelength bifurcation has been identified in an isotropic 
linear system, hexagonal order in the resulting nonlinear structures is actually quite generic HI El . To that 
end, we here consider the model 



2 Model 



dh 
dc 

at 



Ac + BV 2 c + CV 2 h - DV 4 h - EV 4 c 



ac + f3V 2 c + 7 V 2 /i - <5V 4 /i - eV 4 c 



(1) 



(2) 
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describing the (linearized) evolution of perturbations to a steadily-eroding, flat film of a binary material at 
its steady-state concentration. Here h represents the deviation of the surface height profile from that of a 
steadily-eroded flat surface, and c represents the perturbation to the concentration of the atomic species 
that is more easily sputtered by the ion beam. Although we have in mind as an experimental reference the 
irradiation of GaSb, equations of this type arise generically in linear analyses of various kinds of binary 
material irradiation, including the sputtering of a binary material 11201 1281 |2"T1 1221 12311 . sputtering of a pure 
material with impurity co-deposition H24L and ion-assisted growth of a binary material [25]. However, no 
existing works included all of the terms in Eqs. ([l])-(|2]), which represent a superset of existing studies. 

Detailed derivations of each of the terms on the right hand side of Eqns.([l])-([2]) are available in references 
H201 1211 |29H . Here, we give here only a brief, qualitative description of the meanings of each term. 

In Eq.Q describing the height field: (1) the term —Ac (with A > by choice) describes the preferential 
sputtering of species c [20], so that elevated concentration of c leads to enhanced erosion; (2) the term 
BV 2 c describes a net mass flux caused by unequal diffusivities between the components; (3) the term 
CV 2 h measures the dependence of the mass flux on the height gradient, and captures the net effects of 
both curvature dependent sputtering and mass redistribution H301 1311 : (4) the term —DV 4 h (with D > 
by definition) describes many kinds of surface relaxation H301 12611 . and regularizes the height dependence 
in the case that C is negative; (5) the term —EV 4 c would describe a net mass flux from regions of fast 
concentration change to regions of slow concentration change via differing atomic mobilities. 

In Eq.Q describing the concentration field: (1) the term -qc (with a > by definition) describes the 
continual resupply of material into the (c-depleted) irradiated film at the bulk concentration; (2) the term 
/?V 2 c describes simple Fickian diffusion (if positive) or spinodal decomposition (if negative); (3) the term 
7V 2 /i describes species-dependent redistribution H21H : (4) The term —SV 4 h would describe preferential 
migration of species c away from regions of high curvature, due to differing atomic mobilites; (5) the term 
— eV 4 c (with e > by definition) describes an energetic penalty for rapid variations in concentration, and 
serves to regularize the equations in the case that /3 is negative (spinodal decomposition) . 

Comment 1. Equations ([TJl-Q have between them fully ten parameters, most of which are difficult to 
measure experimentally, and nontrivial to measure even with atomistic simulation (see, e.g. H321I331 ]341). 
To slightly simplify our task we follow the reasoning of Bradley and Shipman [12111 and eliminate the highest- 
order cross-dependencies between h and c by setting E = 5 = 0. We also assume here that the diffusivities of 
the two species are equal, which causes B — as seen in H201 1211 . The resulting seven terms still represent 
a considerable number, but we argue strongly that this represents the minimum number of terms required 
to address the question of whether the instability is driven by morphology or composition. In particular, A, 
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D, a, and e are all positive by definition, leaving only the signs of C, (3, and 7 undetermined. We therefore 
anticipate that system behavior will be especially sensitive to the signs of these three parameters. Finally, 
despite the number of parameters, we intentionally avoid the common procedure of reducing their number 
via non-dimensionalzation: despite its appeal, this can obscure implicit sign assumptions on the one hand, 
and physical interpretation on the other. 



3 Stability Analysis 

To investigate the stability of the system, we perform a linear stability analysis, introducing an infinitesimal 
sinusoidal perturbation away from steady state to both the height and concentration fields. Because the 
system is isotropic, we may without loss of generality orient this instability in one co-ordinate direction only; 
we choose the ^-coordinate here, giving 

h 

c 

where k is the wavenumber of the perturbation, a (k) is the dispersion relation for the growth rate of that 
perturbation, and hi and c% are undetermined constants describing the relative phases and magnitudes of 
the height and concentration modulations. We then insert the perturbation ^ into the governing equations 
Q Q in order to find the dispersion relation a (k). 

The value of k that maximizes a (k) - denoted k* - is called the most unstable mode; if a (k*) < 0, then 
the system is stable, whereas if cr(k*) > 0, the system is unstable. In between these regimes is the case 
a (k* ) = 0, which describes a set of points in parameter space called the neutral stability boundary, at which 
the system changes from stable to unstable in a bifurcation. If the bifurcation occurs for k* = 0, it is said to 
be longwave, because the wavelength of the most unstable mode approaches infinity. However, if it occurs 
for some finite value k* ^ 0, the bifurcation is said to he. finite-wavelength, since the most unstable mode has 
a finite wavelength of A = 2n/k*. Generically long wave bifurcations lead to disordered kinetic roughening 
phenomena, whereas finite-wavelength bifurcations lead to ordered patterns H ■ 

Inserting the perturbation Q into the governing equations and ((2), we obtain 



hi 

Cl 



ikx-\-<7t 



(3) 



a hi 



-Aci 



Ck 2 hi - Dk 4 hi 



aci — —aci — f3k 2 ci — ^k^hi 



ek A ci 



(4) 
(5) 
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or in matrix form, 



a + Ck 2 + Dk A A 




hi 


jk 2 a + a + (3k 2 + efc 4 




Cl 



This equation has non-trivial solution for [hi, ci] T only if the determinant of the matrix is zero; which 
condition gives the dispersion relation a (k) - the growth rate a of a Fourier mode of wavenumber k - in the 
quadratic form 

a 2 + T<7 + A = 0, (7) 
where r (k) and A (k) are respectively the trace and determinant of ([6]) when a = 0: 

r{k) = a+(C + ^)fc 2 + (Z) + e)fc 4 (8) 

A (A) = (Ca- A-f)k 2 + (Cf3 + Da)k 4 + (Ce + Dp)k 6 + (De)k 8 . (9) 



For the purpose of identifying stability properties, it is sufficient to consider only the root associated with the 
'+' sign in the quadratic formula which, following common practice, is denoted a + (k): 

v+ (AO = \(- T + Vr 2 - 4A) . (10) 

This gives the growth rate of the faster-growing solution. 

Given the quadratic nature of the equation for a + (k), it is in general not possible to obtain a simple 
closed-form expression for the most unstable mode k*. However, from the form of Equation |10[ we see that 
the system is stable if both t > and a > 0, but unstable if either one is negative. Hence, a bifurcation from 
stability to instability occurs when either 

minr(fc) = (11) 

k 

or 

minA(fc)=0. (12) 

Hence, it is sufficient to examine these two criteria independently. As described elsewhere, both in general 
|[5Tl and in the context of ion irradiation [25], these two possibilities produce different kinds of instability: if 
t (k) > but min/j A (k) < 0, the growth rate is real and the instability is stationary; however, if A (fc) > 
but minfc r (fc) < 0, the growth rate is complex and the instability is oscillatory. 
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3.1 The limit k —> oo: well-posedness 

In part to justify the inclusion of so many terms in our linearized Eqns. ([l]) - ^ , it is instructive to consider the 
limit k — > oo. All well-posed physical problems should exhibit lim / t_ s . 00 a + (k) < 0; i.e., that very-small scale 
perturbations are stable, and consequently, the system will not exhibit any significant variations on scales 
smaller than a certain threshhold. If the limiting value were positive, it would imply the presence of variation 
on infinitely-small scales, violating the continuum hypothesis and suggesting that additional physical terms 
were necessary in the model. 

From Equations (|8])-( 10 1 we quickly find that 



lim a + (k) s» —2ek ; 

k— >oo 

because e > by definition we see that our equations are well-posed. If we were to neglect the term — eV 4 (/>, 
then we would find instead that 

lim cr+ (k) « -2 (D(3) fc 2 ; 

k— >oo 

this is well-posed only if (3 > 0. Hence, the case of phase separation (/3 < 0) is not well-posed without the 
inclusion of the term -eV 4 0. 



3.2 The limit k — > 0: behavior of long waves 

It is also highly useful to examine the limit k — > associated with long-wavelength perturbations. Examining 
Eqns. (|8])-(|1Q|) in the limit as k — >• 0, we find that for small values of k 

lima + (fe)^- (QiC ~ 7A) fc 2 . (13) 

Notably, Eqn.(|13)) depends only on the coefficients of the lowest-order derivatives on h and c, and would 
therefore continue to hold even if no terms had been neglected from Equations Q-Q. We recall that 
A measures the extent to which the component c is preferentially sputtered, and that 7 measures the ex- 
tent to which this same component is preferentially redistributed. If atoms of the preferentially- sputtered 
component are also preferentially redistributed when not sputtered, then A'y > 0, leading to a roughening 
effect, in which c atoms are driven into the valleys (c t ~ 7V 2 h), which then experience enhanced erosion 
(ht ~ —Ac). If this effect is stronger than a smoothing driven by the average rate of mass redistribution 
(h t ~ CV 2 h), then then small values of k have positive growth rates; hence, there is a longwave instability 
leading ultimately to kinetic roughening phenomena. 
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This generic criteria can occur in various physical contexts. For example, in Ref. [24], an increasing 
impurity codeposition rate leads to increased values of 7 relative to C, and the instability criteria Eq. (25) 
there is the condition > aC. Alternatively, in Ref. [12311 . small amounts of the heavier component of a 
binary compound are posited to enhance the sputtering of the lighter, again increasing A'j relative to aC, 
and the instability criteria Eqn.(16) there is again 7 A > aC. Finally, this instability is also documented 
in Ref. [12511 as it applies to the stationary bifurcation in a growing film ( |12| ). Due to its ubiquity and to 
simplify future discussion, we here propose a name for this instability. Because ^A > indicates that 
all displacements - either sputtering or redistribution - occur preferentially for the same component, we 
suggest the name Preferential Displacement Instability (PDI). 

3.3 Finite -Wavelength Requirement and Implications 

Here, however, we are here concerned primarily with the formation of ordered patterns arising when the 
instability has a finite wavelength; the previous discussion serves mainly to reveals that a necessary condition 
for such behavior is the suppression of the PDI instability via 

Ca > A 1 . (14) 

Under the generic finite-wavelegnth restriction ( |14| ), and with the parameters A, a, D, and e all positive 
by definition, a further consideration of the polynomials for r (k) and A (k) reveals that a finite-wavelength 
instability is possible only if one or both of C and j3 are negative. These are notable criteria, because 
C < implies the existence of a morphologically-driven instability even in the absence of concentration 
modulations (ci = 0), whereas /3 < admits the possibility of a compositionally-driven instability even in 
the absence of surface undulations Qi\ =0). In what follows, for clarity and the development of intuition, 
we will discuss the opposing single-mechanism regimes C < 0, (3 > and C > 0, /3 < emphasizing the 
contrast in both instability mechanism and experimental signature. The final regime C < 0, ft < 0, in which 
both mechanisms are operating, is not considered. 

4 Two Distinct Routes to Finite-wavelength bifurcation 

We here present a qualitative discussion of two distinct regimes of finite-wavelength bifurcation. Because 
Eqn. ([9]) is eighth-order, exact bifurcation criteria are rather complicated. Here, we focus instead on differ- 
ences in the signs of key parameters, and the effect on the relative signs of hi and c\ in Eqn. ([6]). 
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(a) 



[1] curvature-dependent sputtering 



A: high-yield 
B: low-yield 




[2] preferential redistribution 




Figure 1 : Representative illustrations of two finite-wavelength instability regimes, (a) In a kinetically-driven 
instability, the surface geometry is fundamentally unstable, while the chemistry follows the interface via the 
preferential redistribution term 7V 2 h in Eqn.Q. To stabilize long waves, it is necessary that the low-yield 
species B be preferentially pushed into the valleys, which leaves the hilltops saturated with the high-yield 
species A. (b) In a chemically-driven instability, the surface is fundamentally stable, following the unstable 
chemistry via the preferential sputtering term — Ac in Eqn.Q. Hence, the hilltops are left saturated with the 
low-yield species B. Mass redistribution acts to stabilize against long waves. 
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4.1 Morphological instability: C < 0, > 

The first finite-wavelength instability regime has recently been discovered by Bradley and Shipman H2 111221 . 
which we briefly summarize for comparative purposes. These authors study the model we have proposed, 
with the exception that they neglect the term — eV 4 c. The lack of this regularizing term, in turn, requires that 
/3 > for well-posedness. Consequently, as described above, instability is only possible if C is sufficiently 
negative. Finally, because these authors sought a finite-wavelength bifurcation, they imposed the finite- 



wavelength criteria (141, which here becomes 7 < —a \C\ /A. 



From these considerations, at the bifurcation point a — 0, the second row of the matrix ([6]) gives 

7 fc 2 /ii + (a + (3k 2 ) 0i = 

Because (a + /3k 2 ) > by assumption, and 7 < is required for a finite-wavelength bifurcation, one 
therefore obtains h\(j)i > 0, so that the concentration variations are in phase with the height variations - 
the preferentially- sputtered species is on the hilltops. The physical interpretation of this scenario is that the 
surface first undergoes undulations due to an erosive instability, after which the preferentially-redistributed 
species is driven into the valleys [12111 . The negative value of 7 implies that the preferentially- sputtered 
species is not preferentially redistributed, and so this species is left on the hills. 

4.2 Chemically-driven instability: C > 0, (3 < 

With the inclusion of the — eV 4 c term, it is no longer necessary to assume that /3 > 0, and here we explore 
the case /3 < 0. This describes a version of the linearized Cahn-Hilliard equation describing phase-separation 
[35], with the addition of a damping term -qc due to material replenishment from the bulk. As motivated 
in the numerical investigation of Le Roy et al. H28L the study of this regime for GaSb is reasonable given the 
presence of the line compound at the 50/50 composition line [36]. In addition, although the low tempera- 
tures typically associated with ion irradiation are far below those needed for thermal segregation, it is known 
that the ion beam induces an effective mobility in the target atoms that is many orders of magnitude higher 
than for an un-irradiated target 026H . Hence, (3 could not only be negative, but of significant magnitude. 
And as a consequence, a finite-wavelength bifurcation no longer requires C < 0. 

Based on these observations, we now consider the case (C > 0, (3 < 0), and re-consider the bifurcation 



criteria ( 11 1-( 12 1. At the bifurcation point a = 0, the first line of the matrix equation ([6]) gives 

(Ck 2 + Dk 4 ) hi + Afa = 0. 
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Because both A > and {Ck 2 + Dk 4 ) > for all values of k, we therefore have h\§\ < 0, so that the 
concentration variations are now out of phase with the height variations - i.e., the preferentially- sputtered 
species is in the valleys. A physical interpretation of this scenario, as suggested in the numerical study of Le 
Roy et al. H28H . is that preferential sputtering of Sb induces the material first to phase-separate into large 
regions of 50/50 GaSb containing small islands of enriched Ga, which then recede more slowly than the bulk 
of the surface. 

4.3 Discussion 

Following Bradley and Harper's seminal demonstration that the erosion process generates a negative contri- 
bution to C H30L it has long been assumed that C < 0, which would yield an unstable surface profile h even 
in the absence of concentration effects. The just discussed erosion-driven instability of Bradley and Shipman 
H211 is consistent with this tradition, and requires concentration effects only to stabilize the long-wavelength 
modes in order to yield a finite-wavelength bifurcation. However, despite its successful reception, we present 
here three reasons for considering alternate mechanisms. 

First, besides atoms sputtered from the surface, the collision cascade associated with ion impact produces 
many atoms that remain on the surface, but are displaced to new locations. The effect of these displacements 
makes a positive contribution to C 113 111 , and very recently, experimental and atomistic studies of pure silicon 
have suggested that in fact this redistribution effect overwhelms the erosive effect, leading to a net value of 
C > H37I 1321 . Although binary systems are chemically different from pure silicon, it is important to note 
that the collision cascade leading to both erosion and redistribution is a purely kinetic process, and hence, 
the kinetic result that redistribution dominates over erosion is likely to persist in binary systems. 

Second, it is frequently observed that irradiated silicon develops a compressive stress H381 1391 1401 1411 
I421I43L and various analytical studies of this effect have recently shown that it, also, makes positive contri- 
butions to C H44I II PL despite initial speculation to the contrary 0E1I- Again, although the analogy with 
binary systems is not ensured, stress generation is likely to persist in binary systems as long as it is caused by 
the kinetic collision cascade. (However, see Ref. [16], in which the chemical reaction of Mo and Si to form 
MoSi2 leads to a tensile stress hypothesized to produce an instability.) 

Third, the specific requirement on concentration effects from the BS theory is that 7 < —\C\a/A, so 
that the species more likely to be sputtered, is less likely to be redistributed if left on the surface. However, 
erosion and redistribution are merely different outcomes of the same kinetic process - a displaced atom 
either has enough momentum away from the surface to escape, or it does not. From this perspective, it is 
not immediately intuitive that a single species would be both less likely to be displaced far from its initial 
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position, but simultaneously more likely to escape the target. 

The significant restrictions on the unknown parameters C and 7 are an important additional motivation 
to investigate the plausibility of phase separation in driving ordered patterns. As we have seen, if j3 is 
sufficiently negative to drive a bifurcation via phase separation, then there is much more flexibility in the 
values of C and 7. First, C could be positive, allowing ordered patterns even if the destabilizing erosion 
process is not strong enough to dominate other, competing stabilizing processes. And second, a positive 
C would also serve to stabilize long-wave disturbances via Eq.(|14P, allowing the possibility that the same 
component is both preferentially sputtered and redistributed. 

5 Conclusions 

The primary contribution of this work is the introduction of a generalized linear framework in which both 
kinetically- and chemically-driven instabilities can be compared analytically. Our main finding is that phase 
separation admits a second finite-wavelength instability, which exhibits the opposite experimental signature 
of the BS instability: the erosion-driven BS instability requiring C < and 7 < exhibits concentration 
modulations of the preferentially-sputtered species that are in phase with the height modulations, whereas 
the chemistry-driven instability driven by (3 < exhibits out-of-phase modulations. Unfortunately, although 
both mechanisms have been plausibly motivated already in the literature H21U28L better data are urgently 
needed to choose between them, and the main utility of our study is to suggest two specific routes to resolving 
the question of which mechanism is dominant. 

Regarding parameter values, existing studies on pure Si H32[ |34B have been successful in comparing the 
relative contributions of erosion and net redistribution to the value of C. Such studies could be general- 
ized for the two-component GaSb system without much alteration. In addition, they could simultaneously 
investigate the value of 7, which describes preferential redistribution. If both parameters were found to be 
negative for GaSb, it would bolster the BS theory, while if either were found to be positive, it would instead 
support phase separation. 

Experimentally, it would be invaluable to determine the relative phases of Ga and Sb at the early stages of 
pattern formation. Notably, disagreement currently exists within the literature. As mentioned above, using 
ex-situ electron energy loss spectroscopy (EELS) applied to a single nanostructure, Le Roy et al. observe 
very well-defined Ga-rich caps at the top of the structure H27H . However, using low-energy ion scattering 
spectroscopy (LEISS) applied to a large array of nanostructures tilted at various angles, El-Atwani et al. 
come to the opposite conclusion, suggesting that the peaks of such structures may instead by Sb-rich [14511 . In 
addition to these conflicting observations, both studies focus on late times, after a well-defined pattern has 
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already developed, in contrast to the early times associated with linear analysis. 
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